Emergence of robust gaps in 2D antiferromagnets via additional spin- 1/2 probes 
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We study the capacity of antiferromagnetic lattices of varying geometries to entangle two addi- 
tional spin-1/2 probes. Analytical modeling of the Quantum Monte Carlo data shows the appearance 
of a robust gap, allowing a description of entanglement in terms of probe-only states, even in cases 
where the coupling to the probes is larger than the gap of the spin lattice and cannot be treated per- 
turbatively. We find a considerable enhancement of the temperature at which probe entanglement 
disappears, as we vary the geometry of the bus and the coupling to the probes. 

PACS numbers: 



Solid state systems have been exploited to integrate 
quantum information (QI) tasks and accomplish quan- 
tum computation in a single processing core, but many 
questions regarding their robustness against temperature 
and decoherence remain open. Among the requirements 
to achieve quantum computation, the ability to gener- 
ate rapid elementary gates between well-characterized 
qubits is central [H . In view of the technical difficulties 
of switching on direct interactions between qubits, vari- 
ous proposals have been put forward to use a quantum 
sub-system, usually denominated as 6ms, to mediate the 
fundamental universal gates. It can be realized, for in- 
stance, by the phononic mode of cold ions, in the famous 
ion-trap quantum computer 2], or by the magnetic de- 
grees of freedom of a quantum spin chain (see (3| and ref- 
erences therein). A considerable body of work has been 
devoted to chains of spins experiencing nearest-neighbor 
interactions, since they can be used as models for univer- 
sal quantum computation, meeting the aforementioned 
requirements 0]. Spin chains have been shown to be 
extremely versatile; for instance, they allow the reliable 
transfer of the state of a single qubit , and their ground- 
states are able to mediate an effective long-distance inter- 
action ultimately entangling distant spin probes — the 
so-called long-distance entanglement (LDE) 0] . 

Promising advances in the engineering of atomic struc- 
tures and optical lattices, where finite spin systems are 
effectively realized in the laboratory, encourage the con- 
sideration of more general possibilities. We consider anti- 
ferromagnetic (AF) spin systems, ranging from ID chains 
to square lattices, and demonstrate the opening of robust 
gaps via local interactions with two additional spin-1/2 
particles. Robust gaps implies negligible thermal occu- 
pation of excited bus states in the entire range of temper- 
atures in which the probes are entangled; this allows for 
QI capabilities {e.g. LDE) at higher temperatures than 
previously considered possible. 

Models of spins — We start by making precise the no- 
tion of "robust gap"; by this we mean that the bus plus 
probes system has a low energy manifold isomorphic to 
the probe space, well separated in energy from the re- 
maining energy spectrum (see Figures [T] and [2|) . Should 
this be so, and for the case of rotational invariant cou- 



Figure 1: The probes (blue) interact locally with the lattice 
through sites A and B. Will the bus create an effective inter- 
action capable of entangling the probes at finite temperature? 



plings, symmetry alone allows for an analytical descrip- 
tion of the temperature dependence of the entanglement 
between probes. Surprisingly, we found that in all the 
systems we studied by means of Quantum Monte Carlo 
(QMC) simulations, entanglement follows this descrip- 
tion very accurately. This shows that, even in systems 
where the bus has a very small gap (approaching zero 
in the thermodynamic limit), the actual coupling to the 
probes can ensure the emergence of a robust gap, an ideal 
situation for quantum computation. 

We consider two qubit probes that interact locally with 
a generic finite spin lattice (the bus). Setting h = 1, the 
Hamiltonian of the full system (bus -I- probes) reads. 



H = Hhus + aJ {Ta ■ Sa + U- Sb) ■ 



(1) 



where Hjyus is the Hamiltonian of the bus, a J is the cou- 
pling strength between the probe qubits and the lattice, 
and J a characteristic energy scale of the bus. Eq. ([1]) 
describes an isotropic interaction between the probes - 
with Pauli operators, T4(b) - and the bus spins at spe- 
cific lattice sites (5^(5); A,B e C). We speak about 
LDE at temperature /csT = whenever the partial 
state of the probes is entangled for distances of the or- 
der of the system size, (Iab ~ 0{L). Genuine quantum 
correlations among qubits living in the bulk system usu- 
ally decay rapidly [8(], and, only recently, one-dimensional 
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Figure 2: Schematic of the opening of a robust gap (A(q) ^ 
fcflT) by two spin- 1/2 probes that couple locally to the many- 
body system (bus) with arbitrary strength aj. This is ob- 
served in all the AF lattices considered in this article. If the 
singlet is localized near the additional probes, they will be 
highly entangled even at large distances - the so-called LDE. 
The singlet has a protection gap (triplet-singlet energy sepa- 
ration, J) which is enhanced with the dimensionality of the 
bus. 



spin systems were found to support robust LDE [1, 0| , as 
opposed to gapless bosonic systems |9|] . 

We focus on SU (2) symmetric Heisenberg interactions, 
which not only allow universal quantum computation Q , 
but are also commonly realized in nature (e.^. in the 
parent compounds of copper-oxide high-temperature su- 
perconductors, such as the undoped insulator La^CuO/^ 
[lOj : in electronically coupled quasi- ID chains such as 
CuGeOs ll|; in the Mott insulating one-dimensional 
perovskite, KCuF^ [13] ; in linear chains of 10 man- 



ganese atoms in engineered structures [13|). 

The partial state pat is completely fixed by symmetry 
and depends on a single physical parameter driving the 
bus capacity to entangle the probes, 



Pab = CXp {-PJabTa ' n) . 



(2) 



The effective coupling changes with temperature, Jab = 
Jab{P), since tracing out the degrees of freedom in C in- 
troduces mixedness. A QMC calculation of the correla- 
tion between probes allows easy extraction of Jab{fi)- As 
it stands, Eq. ^ would mark no real advance, if not for 
the fact that symmetry allows a complete specification 
of the temperature dependence of the effective coupling 
Jab = Jab{P), provided the condition of a robust gap is 
verified; Jab can be expressed as function of three tem- 
perature independent parameters {J , and 77) in the 
entire parameter region where the probe entanglement is 
not too small: 



Ja6(/?) = ^ln 



3($ - ?7) (4 - 3$ - r/)e' 
4-$-|-77-|-($ + ?//3)e'3-^ 



(3) 



The canonical coupling is the gap between the low- 
est singlet and triplet (the protection gap), and sets 
the temperature scale for the disappearance of entangle- 
ment. The canonical corrections, $ and 77, determine the 
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Figure 3: Black triangles: effective coupling. Jab, at a dis- 
tance I — 20, as a function of ric, number of transverse chains 
{ksT — 2 X 10~^J ); blue dots and red squares: critical tem- 
perature above which LDE vanishes; inset: canonical coupling 
(i.e. the probes singlet protection gap of Fig. [2]), in the same 
units used to represent Jab- The error bars from QMC cannot 
be seen as they are typically below 1%; a = 0.05 in all plots. 



amount of entanglement at zero temperature (see Ap- 
pendix). Although these parameters cannot be analyt- 
ically computed in general, they are well-approximated 
by their values in perturbation theory for a sufficiently 
small. For weak coupling, it turns out that since 
the first non-zero contribution to rj is of order a^, whereas 
$ ~ 0{a^). The formula ^ is exact and will be valid 
whenever the bus-probe system has a robust gap. 

Numerical results: Effective interaction enhancement 
in 2D — We now present the main results from the QMC 
simulations on a family of AF lattices and compare the 
results with our theoretical prediction [Eq. ^] . Our sys- 
tems consist of 2-dimensional (2D) finite lattice £, with 
N = I X Uc spins-1/2 and two extra probes, where / is the 
number of longitudinal sites and stands for the num- 
ber of coupled chains, varying from — I to — I. The 
Hamiltonian of the lattice is Hq = J^i^ij^ Si ■ Sj, with 
J > 0. The qubit probes interact with the spins at the 
boundary of the most central chain (see Fig. [1]) through 
an isotropic interaction [Eq. We expect a significant 
change in the LDE from the common one-dimensional 
scenario [01, as the physics of a 2D bus is very distinct. 
In particular, the 2-leg ladder chain has an Haldane gap 
which should play against a large Jab since very mas- 
sive excitations, A ~ 0.504J, make the correlations die 
particularly fast 
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We briefly present the results for small probe-bulk cou- 
pling, a ~ 0.05, before venturing away from perturba- 
tion theory. The results of Fig. [3] and Fig. [4] show a 
clear enhancement of the ability of the antiferromagnet 
to generate long-range effective interactions among dis- 
tant probes as one reaches the square-lattice. Also, a wig- 
gly behavior up to nc = 4 is evident (the large Haldane 
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Figure 4: The points in the plot show Jab as function of the 
temperature from the QMC simulations for a — 0.05. The 
hnes stand for the fit with the expression given in Eq. |[3]). 
All curves saturate for high-temperatures [Eq. Q] and there- 
fore entanglement vanishes much before our method becomes 
inaccurate. The agreement between the QMC data is excel- 
lent resulting in a average deviation of ~ 0.1 — 1% depending 
on the lattice. 



gap for He = 2 and 4 implies a fast decay of correlations 
within the bulk spins leading to a negligible effective cou- 
pling). Fig. [3] shows a curious transition for ric > 4: the 
increase of the protection gap (and also Jab) becomes 
smooth and the Haldane finite-size gap, very strong for 
TT-c = {2, 4}, gets suppressed — the 2D physics is reached 
monotonously as the bus gap disappears. 

We expect the ground state of 2D antiferromagnets 
to reduce substantially the LDE due to the symmetry 
breaking at T = 0, for large lattices; the finite stag- 
gered magnetization should reduce the amount of gen- 
uine quantum correlations shared by the probes. This is 
borne out by the results of the QMC simulations, shown 
in Fig. 21 where Jab is found to decrease at low tem- 
peratures, when the number of chains increase. Never- 
theless, at higher temperatures, the opposite occurs. Jab 
increases with ric] this reflects the increase of the protec- 
tion gap, J^. Having shown the QMC results for 20 spin 
lattices, we now compare them with Eq. ([3]) for several 
temperatures: Fig. H] shows a perfect fit to the QMC 
data. For sake of clarity, we have presented the agree- 
ment just for 4 lattices although all them show the same 
degree of accuracy. The observed linear dependence of 
Jab with the temperature for T ^ is easily understood: 
a zero temperature (finite) entanglement below the max- 
imum value of 1 requires Jab/T constant (see Fig. 
[5]). This constant can be derived from Eq. ([3]), yielding, 
4/3Jah In [(4 - 3$ - ?7)/($ + 77/3)] ; thus, the canon- 

ical corrections (<i> and 77) determine the low-temperature 
physics of the probes. For high-temperatures the effective 
couphng satures (see Fig. 2]) to a constant value, ^7(1 — $), 
when r] is negligible {i.e. not far away from the perturba- 
tion limit, see Eq. ^ and comments therein), suffering 
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Figure 5: Concurrence, an entanglement monotone for qubits, 
as function of temperature for representive lattices and dif- 
ferent couplings. Lattices supporting more entanglement at 
strictly T — have worse performance at higher tempera- 
tures. 



a slightly increase otherwise. 
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(4) 



Figures [3] and |4] deal with relatively small probe-bus 
coupling, but the results presented so far are more gen- 
eral. For instance, choosing a sufficiently large a to 
strongly suppress the zero temperature entanglement, 
due to partial frustration among the neighborhood of the 
bulk spins connected with the probes, we again find an 
excellent agreement with Eq. ([3]); the measured concur- 
rence is perfectly fitted with an expression derived from 
Eq. ([31) (see Fig. [3). This clearly shows that in all our 
measured systems, the condition of a robust gap is ver- 
ified. This is surprising, particularly in the case of the 
(finite) square lattice, which has a gap much smaller than 
aJ; one would not expect, in this situation, the appear- 
ance of a well protected singlet. The emergence of the 
robust gap has to be attributed to the coupling of the 
probes: the lowest singlet and triplet are pulled down 
from the rest of the spectrum, allowing a complete de- 
scription of entanglement only in terms of these two en- 
ergy levels. On the other hand, whereas this stronger 
coupling to the bus reduces the zero temperature entan- 
glement, it also allows a larger split between the singlet 
and triplet, leading to entangled probes at much higher 
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temperatures. Typically, exchange interactions in anti- 
ferromagnets can be of the order of 0.1 eV, resulting in an 
effective coupling of the order of 0.3 meV for the square 
lattice {I — 20) at temperature ~ 12K and a = 0.2. 
This is to be compared with the value of 0.01 — 0.1 meV 
achievable in quantum dot spins [3] (see comment [13] )■ 
Regarding LDE (see also Fig. [5]) we see that the criti- 
cal temperature (above which the correlations shared by 
the probes are completely classical) can be increased by 
a factor of 20 from weakly coupled spin chain {a = 0.05) 
to an intermediate coupled (a = 0.2) square lattice, en- 
tanglement surviving up to ksT ~ 1.2 x 10^^ J. 

Conclusions — We observed an enhancement of LDE 
by varying the geometry of the magnetic spin systems 
serving as a quantum bus. The analysis based on the 
canonical transformation formalism agrees very well with 
the results from QMC simulations, demonstrating the 
emergence of robust gaps in AF lattices via a coupling 
of intermediate strength with external probes. We raise 
the possibility of entangling distant spin probes at tem- 
peratures as high as T ~ 1.2 x 10~^J//cb, where J is the 
nearest neighbor exchange constant of the bus. We hope 
this work stimulates research in the area of "all-in-one" 
solid-state based quantum communication and computa- 
tion. 

Appendix: Analytic treatment and the canonical pa- 
rameters — We take the bus to be a spin 1/2 spin system 
with rotationally invariant couplings and a non degener- 
ate singlet ground state |0o); the probes couple to the 
bus also through isotropic couphngs, Eq.([T]). If adiabatic 
continuity holds ^16*1, we have a one-to-one map of the 
eigenstates of the full Hamiltonian to those of a = (un- 
coupled probes) ; this map defines a canonical transforma- 
tion, \(j}m) ® IXcr) = e*'^IV'm,o-)j where |0,„) is a bus-only 
eigenstate, |xcr) a probe state, and IV'm.cr) an eigenstatc 
for finite a. Under this transformation, the Hamilto- 
nian must separate into a bus term and a probe term, 
since the corresponding eigenstates are product states. 
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of this canonical transformation is generally achieved in a 
perturbation expansion in a [l^ , but we choose to keep 
it unspecified for the moment. If we use the canonical 
transformation in the expressions of thermal averages, 
we can easily rewrite them in terms of Tis', on the other 
hand, the assumption of a robust gap, i.e., the lowest 
lying states of the full system, which map to the product 
of 100 ) and the probe singlet or triplet, are well separated 
from remaining states by a gap A(a) ^ fcsT, allows the 
neglect of thermal occupation of states which map to ex- 
cited states of the bus; one obtains 



{Ta ■ Tb) 



Tr, [e-P^^Y^rnAraTa-TbAl] 

Trp [e~P^^] 



where Trp(. . . ) is a trace over probe states, and Am — 
{(j^oW^ \'t'm) is a probe-only operator. Rotational symme- 



try implies that A^Ta ■ TbAl^ must be a scalar oper- 
ator [Q]; therefore, J^m. ^mTa ■ nA]^ = jy-^ (1 - $) Xa ■ n, 
where, by construction, r] and $ are temperature inde- 
pendent. This obviously entails (Xa -Tb) = ry-l- (1 — <I>) (tq ■ 
'ni)can, where (. . . )can IS an average with the Hamiltonian 
Tip, which has the form H. = \jTa ■ Tb, implied by sym- 
metry. By expressing (xa • Tb) in terms of Jab, using 
Eq.(l2]), and (tq • 'n,)can in terms of J, one easily derives 
Eq. Q. 

The actual calculation of these parameters, J , $ and 
?7, can only be done in perturbation theory, following the 
approach of Schrieffer and Wolff [19]. The calculation 
of the effective coupling J has been discussed in [7|; to 
calculate the canonical corrections 77 and $ we start from 



boW^Ta-Tbe-'''\^o) 



Ta ■ Tb+ liCpol 



S,Ta ■ Tb 



S, Ta ■ Tb 



and use the generator prescribed by Hp-b + iS,Hb =0 
[19I, to obtain 



rj = + O{a 



m>0 
4^ 



1 



{Em — Eq) 



\{MiSl-SI,)\cPm)\}{5) 
(6) 



The fact that 77 is of higher order than $ means that, for 
weak coupling, we have a simpler result than would be 
implied by rotational symmetry alone, namely, (xa - Tb) = 
(1 - $)(ra • Tb)can, and Jab giveu by Eq.© with 77 = 0. 
In fact, for a < 0.1, a careful inspection shows that the 
fits we present are virtually indistinguishable from the 
fits with rj = 0. 
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